Current fluctuations in composite conductors: Beyond the second cumulant 
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Employing the non-linear cr-model we analyze current fluctuations in coherent composite con- 
ductors which contain a diffusive eiement in-between two tunnel barriers. For such systems we 
explicitly evafuate the frequency-dependent third current cumulant which also determines the lead- 
ing Coulomb interaction correction to shot noise. Our predictions can be directly tested in future 
experiments. 



I. INTRODUCTION 

Investigations of current fluctuations in mesoscopic 
conductors enable one to extract important information 
about their microscopic properties, such as, e.g., trans- 
missions Tk of their conducting channels. A great deal of 
information is provided by the analysis of shot noise^. 
For instance, experimental confirmation of shot noise 
suppression in diffusive conductors (as compared to its 
Poisson value) by the Fano factor (3 — 1/3 served as a 
strong argument supporting the validity of the universal 
transmission distribution P{Tk) oc 1/Tk\/1 — Tk in such 
conductors, see Refs. [lll^for more details on this subject. 

It was also realized, that transmission distributions in 
composite conductors, such as, e.g., a tunnel junction 
attached to a diffusive wire can behave in a rather non- 
trivial way. For instance, it was observed experimen- 
tally that the low energy Andreev conductance of hybrid 
normal-superconducting structures may coincide with its 
conductance in the normal state. This zero-bias anomaly 
phenomenon^ was explained^i^ by finding out that the 
transparency distribution for a composite system reduces 
to that of a diffusive conductor provided the resistance of 
the latter exceeds that of the tunnel barrier. This phe- 
nomenon can be interpreted as disorder-induced opening 
of tunneling channels. 

Yet another non-trivial property of composite conduc- 
tors involving many coherent scatterers is that shot noise 
in such systems turns out to be identical to that in a 
diffusive element even if none of these scatterers is diffu- 
sive. These property was initially established for chains 
of tunnel barriersS, and later demonstrated also for other 
situations^i^. Most generally, one can consider a chain of 
N different coherent conductors with arbitrary transmis- 

(n) 

sion distributions T^, of their conducting channels and 
define their resistances i?„ and Fano factors /3„ as 
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with e standing for the charge of the electron. As usu- 
ally, shot noise in such composite conductor with resis- 
tance R = ^n=i '"^^ described by the correlation 
function 



where SI{t) = I(t) — I accounts for current fluctuations 
around its average value I = V/R and V is externally ap- 
plied voltage bias. Here and below averaging (...) should 
in generally be understood as the expectation value for 
the corresponding current operators. For the symmetric 
version of the correlator ^ in the zero frequency limit 
and at zero temperature one has 52(0) = |e/|/3, where 
the Fano factor (3 of this general composite conductor 
reads^ 
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This formula clearly demonstrates that in the limit of 
sufficiently large N one has (3 — >■ 1/3, i.e. a long chain of 
arbitrary - not necessarily diffusive - conductors univer- 
sally behaves as a diffusive conductor. 

Along with shot noise one can also investigate higher 
current cumulants. According to full counting statistics 
(FCS) theory^ these cumulants are in general expressed 
not only in terms of the Fano factor (3 but also con- 
tain other combinations of channel transmissions. For 
instance, the third current cumulant 

cS3(wi,W2) - J dndT2e''^i^i+^"=^^ 

x{SI{t)6I{t-n)SI{t^T2)) (4) 



after its proper symmetrization^*^ can be cast to the 
formiajJ- 



53(wi,W2) = (/3-27f^)e2j, 



where 



7 ^ 
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52 (w) = / dre"^^ {5I{t)5I{t - r)) 
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and F is the function of bias voltage V , temperature 
T and frequencies aji.2 to be defined below. Thus, by 
studying the cumulant (|4]) one can recover not only the 
Fano factor (3 but also the combination ([5]). Experimen- 
tal investigations of the third current cumulant have al- 
ready startedi2,j thereby making 7 (along with (3) a di- 
rectly measurable quantity. Furthermore, in the presence 
of electron-electron interactions the combination ^ also 
affects the second current cumulant^. Hence, 7 can in 
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principle be extracted from the noise measurements as 
well. 

In the light of all these developments it would be useful 
to formulate a general approach enabling one to recover 
the parameter 7 ^ for composite conductors. In the 
present paper we will address a specific situation that is 
frequently encountered in modern experiments embrac- 
ing a variety of meso- and nanostructures: An arbitrary 
coherent diffusive conductor is connected to an exter- 
nal measuring scheme (leads) via tunneling interfaces. 
In other words, we will consider a composite conductor 
consisting of two tunnel barriers connected by a diffusive 
element. Employing the non-linear a- model technique we 
will develop a method enabling one to evaluate the FCS 
generating function for this system. We will then derive 
the parameter 7 in terms of both tunneling and diffusive 
resistances involved in our problem. The corresponding 
expression can be directly used in future experiments in- 
vestigating current fluctuations in such composite con- 
ductors. Conversely, information about three first cur- 
rent cumulants allows to fully determine the resistances 
of the three elements of the system under consideration. 

Our paper is organized as follows. In Sec. II we work 
out a non-linear cr-model approach and establish general 
expressions for the parameters /? and 7 of our composite 
conductor. Our results and their possible extensions are 
then discussed in Sec. III. Some technical details of our 
calculation are relegated to Appendix. 



II. EFFECTIVE ACTION AND FCS 

We will consider a system which consists of a normal 
diffusive conductor with resistance Rd and length L at- 
tached to two big normal reservoirs via tunnel barriers 
(located at x = and x = L) with resistances Ri and R2 
and cross-sections Fi and This system is schemati- 
cally depicted in Fig. 1. We will assume that the whole 
system remains coherent, i.e. L is shorter than both 
phase and energy relaxation lengths. We will also as- 
sume that the Thouless energy (inverse dwell time) of a 
composite conductor in-between two reservoirs is larger 
than any other energy scale in our problem. 

In order to proceed below we will employ the Keldysh 
non-linear cr-model approach-^. The kernel of the evo- 
lution operator defined on the Keldysh contour can be 
expressed in terms of the path integral 



j VgVAVV eM^Si^.\ds + iS), 
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where A and V are respectively the vector and scalar 
potential fields and g{R,ti,t2) are 2x2 matrices which 
depend on one coordinate and two time variables and 
obey the normalization condition 



j dt'g{R,ti,t')g{R,t',t2) 



S{h-t2). 
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FIG. 1: Diffusive coherent conductor with resistance Rd be- 
tween two normal reservoirs attached via tunnel barriers with 
resistances Ri and R2. The barriers are located at a; = and 
X — L. 



Here S'fioids[^, V^] is the action for the electromagnetic 
fields and 



S — Sd Si + S2 



(9) 



defines the action for our composite conductor. Specifi- 
cally, 



S 



iirNo 



D 



Tr 



D {dgf + Ai {idt - eV^ g 



(10) 



is the non-linear cr-model effective action which accounts 
for a diffusive conductor and the terms 



^1 



TTCTi 



Trri 513(0), ^2 



7rcr2 



Trr,g(L)g2 (U) 



describe tunnel barriers. Here and below cri_2 = 
1/(^1. 2F1, 2) and gi^2 represent the quasiclassical energy- 
integrated Green-Keldysh matrices of the reservoirs to 
be specified below. The expressions pT|) follow di- 
rectly from Kupriyanov-Lukichev boundary conditions^—. 
These boundary conditions are sufficient to employ here 
since our further analysis will be restricted to the lowest 
order in barrier transmissions. 

Here and below the product of matrices g should be un- 
derstood as a convolution, cf. Eq. ([8]). The trace is taken 
over both space and time variables and it is accompanied 
by a summation over matrix indices. The diffusion co- 
efficient and the density of states per spin at the Fermi 
surface of a diffusive conductor are denoted by D and A^o 
respectively. Spatial integrations in the terms 6*1,2 (fTTj) 
are restricted to the respective barrier interfaces. 

The 2x2 matrix V term depends on the Hubbard- 
Stratonovich fields and defined respectively on 
the forward and backward branches of the Keldysh con- 
tour. As usually, these fiuctuating fields emerge after the 
standard decoupling procedure in the Coulomb term in 
the Hamiltonian. Here we define 



V{R,t) 



2 ^ 



(12) 



where y+ = iy^ + V^)/2,V- ^V^ -V^ are respec- 
tively the "classical" and "quantum" fluctuating fields. 
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Finally, we specify the operator 



dg = Vg - 



A 9 



(13) 



Wg{t, t') ~ i-A{t)g{t, t') + g{t, t')i-A{t'). 
c c 



In Eq. (|13p we suppressed spatial arguments for the sake 
of brevity, introduced the speed of light c and defined the 
matrix 



1 A 



(14) 



is the Pauli matrix. The transformation (jl7l) yields g2 in 
the form 



he 



(19) 



cos Lp -\- i sin ipN{e) 
-i simp 



2sm^ fiV(e) 



isin(^ + 2cos^ ^N{e) 
— cos Lp — i sin ipN(e) 



What remains is to evaluate the path integral ([7]) 
within the saddle point approximation. The saddle point 
paths obey the well known Usadel equatio n^^i^^ . In the 
stationary limit to be considered below this equation 
reads 



which - analogously to Eq. ([T2|) - is composed of the 
classical = {A^ + A^)/2 and quantum A = A^ — 
A^ components of the vector potential. 

Now let us employ the above non-linear cr-model for- 
malism in order to construct the PCS generating func- 
tion. In what follows we will assume that an external dc 
voltage V is applied to the system. This voltage is ac- 
counted for by a simple shift of the classical component 
V"+ V + . At this stage we also neglect the effect 
of electron-electron interactions by setting the quantum 
component of the scalar potential equal to zero V~ — >■ 0. 
At the same time we will still keep the quantum com- 
ponent of the vector potential A~ which is needed for 
further evaluation of the current cumulants. Note that 
the frequency dependence of the third current cumulant 
in Eq. ([5]) is fully accounted for by the function F eval- 
uated in^*', while the parameters /3 and 7 are frequency- 
independent. Hence, in order to determine these param- 
eters it suffices to perform our calculation in the zero 
frequency limit. In this case, similarly to Ref. [l^ it is 
convenient to employ the pure gauge 



(15) 



We will see that x is going to play the role of the counting 
field in the FCS generating function for our problem. 

Now let us specify the Green-Keldysh matrices of the 
reservoirs. Setting the scalar potential of the left reser- 
voir equal to V we have in the frequency representation 



5le 



1 2N{e-eV) 
-1 



A^(e) = tanh(e/2r). (16) 



Accordingly, the scalar potential of the right reservoir 
equals to zero. The counting field x can be gauged away 
from the action (|10p at the expense of the following trans- 
formation of the Green-Keldysh matrix of the right reser- 
voir: 



32 new = exp 
where (p = exi c and 



yCTj. 1 52 old exp I yCTa 



1 

1 



(17) 



(18) 



D[d,g,{x)[d,g,{x)]] = [-ie.g.i^)]- 



(20) 



Eq. (I^Hl) should be supplemented by Kupriyanov- 
Lukichev boundary conditionsi^ involving the Green- 
Keldysh matrices of the reservoirs ([TS]) and Once 
the solution of Eq. (|20p is established, the generating 
function for current cumulants iS{x) is immediately ob- 
tained by substituting this solution back into the action 

®- 

Let us define the matrix 



J(e) = gedxge 



(21) 



which anticommutes with g^^ and does not depend on the 
coordinate x. The latter property just reflects the spec- 
tral current conservation. With the aid of this matrix the 
solution of Eq. (^01 inside the diffusive conductor can be 
expressed in the form 



geix) = ge(0)exp (^J(e)i 



(22) 



In addition, Kupriyanov-Lukichev boundary conditions 
yield 



J(e) 



2<JD 
2CTD 



[3ie,.9.(0)], (23) 
[ge(0)exp(J(e)L),g2e], (24) 



where ao — 2e^DNo is the Drude conductivity. 

Let us fix some parametrization of ^£(0) satisfying the 
normalization condition g^{0) = 1. Then with the aid of 
Eq. ([23|) we evaluate the matrix J(e). Substituting the 
resulting expression into Eq. (j24p we determine the pa- 
rameters employed in our parametrization for ge(0). Fur- 
ther details of this calculation are presented in Appendix 
A. Here we evaluate the effective action up to the third 
order in x- Making use of the identity Tr {dxgf' — — Tr 
and keeping the contribution from the interface terms we 
eventually arrive at the generating function 



■Q( \ 

c 



elx^ 
2c2 



P coth ^ 



2(1-/3) 



ie^lx^ 
6c^ 



{(i-2^F{v)) + 0{x^), 



(25) 
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where 



1 , 2{Rl + Rl) 



3 3 i?3 

2 6(i?f + i?l) 4(i?f + i?^)2 



7= 77 + 7 



15 5 i?5 3 i?6 



(26) 
(27) 



Here we defined v = eV/T, the total resistance R 
Ri+Rd + i?2 andii 



F{v) = 1 + 3 



1 — (smhv/v) 
cosh V — I 



(28) 



The generating function (l^5|) - (l?fl) represents the central 
result of this paper. Making use of the correspondence 
I — >■ —ic5/5x this result allows to recover the first three 
current cumulants for our composite conductor. It is im- 
portant to emphasize that the expressions for these cu- 
mulants - including the third onei^ - can be fully de- 
scribed not only in the zero frequency limit, but also at 
finite frequencies. This and some other applications of 
the above results will be briefly discussed in the next 
section. 



III. DISCUSSION AND OUTLOOK 

It is satisfactory to observe that our result for [3 ((26)) 
is fully consistent with the general formula ([3|) derived 
within a totally different framework. We also observe 
that, as compared to /3, the expression for 7 ([27]) is deter- 
mined by a substantially more complicated combination 
of resistances i?2 and Ru. Even more complicated 
expression for 7 can be expected in a general case of N 
arbitrary conductors. 

Let us specify our results for /3 and 7 in several impor- 
tant limits. Provided the resistance of one tunnel bar- 
rier dominates i?i > Rd,R2, Eqs. dMl) and d^T]) yield 
/3 « 1 and 7 « 0, i.e. in this limit current fluctuations 
are identical to those in a tunnel barrier. If, however, 
the diffusive resistance prevails Rd 3> i?i , i?2 , we have 
/3 w 1/3 and 7 w 2/15. These values are characteristic 
for any diffusive system. In the case i?i ~ R2 ^ Rd 
from Eqs. dMl and (HZ]) we obtain ;3 « 1/2 and 7 « 1/8. 
The latter value of 7 turns out to already be rather close 
to the diffusive one 7 = 2/15 rapidly approaching it with 
increasing Rn- E.g., for Ri ^ R2 ~ Rd the Fano factor 
/3 still remains bigger than 1/3 by 15% while 7 differs 
from 2/15 by less than 2%. 

Turning to the frequency dependence of the third cur- 
rent cumulant and making use of the results^*^ we estab- 
lish the function F in Eq. ([S]) in the following general 
form 



F{v, wi, W2) = Z{v, ~wi,wi + W2) 

+ Z{v, -W2,Wi + W2) + Z{v, 'Wi,W2), 

where wi^2 = '^i,2/2T and the function Z reads 
sinh(u/2) 



Z{v,x,y) 



2v sinh x sinh y 



(29) 



(30) 



x + y-{v/2) 



x + y + {v/2) 



smh[x + y — {v/2)] smh[x + y + {v / 2)] J 

In the zero frequency limit W1.2 — >■ the function F 
reduces to that defined in Eq. (^51) . The above equa- 
tions together with Eqs. and ([271) fully describe 
the third current cumulant in our composite conductor 
provided its size is restricted in such a way that its Thou- 
less energy (inverse dwell time) remains the highest en- 
ergy parameter in the problem. Note that the opposite 
physical limit of very lone dwell times was studied for 
a chaotic cavity in Ref. il9l The frequency dependence 
of 53 established in that limitiS differs from ours, how- 
ever in the limit of zero frequencies the result o^^'^*^ are 
consistent with ours (taken at T — ^ 0) provided we set 
Rd ^ Q in Eqs. (gS]) and (gT]). We also note that non- 
symmetric third cumulants of the current can also be 
analyzed^i. The corresponding results for 53 being in 
general different from Eq. ([5]) also contain the parameter 
7 (O evaluated here. 

Let us now take into account electron-electron inter- 
actions. It was argued", that such interactions may af- 
fect shot noise. In particular, for R <^ Rq and T,\uj\ <C 
|ey| < l/RC we obtaini^ 



S2{uj) = \el\p- 



2{(3-2-t)\eV\ 

Rn 



In 



1 



eV\RC' 



(31) 



where Rq = h/e^ and C are respectively the quan- 
tum resistance unit and the effective sample capacitance. 
This result also follows directly from the renormalization 
group analysis22i2^ stating that electron-electron interac- 
tions effectively yield energy dependent renormalization 
of channel transmissions. Eq. ([31]) demonstrates that 
the magnitude of the leading interaction correction to 
the second current cumulant 52 is governed by the com- 
bination j3 — 27. This combination - depending on the 
system - can take either positive or negative values. Ac- 
cordingly, electron-electron interactions can either sup- 
press or enhance shot noise. In the situation considered 
here the parameter /? — 27 is always positive ranging from 
1/15 to 1 depending on the resistance values. Hence, in 
our case Coulomb interaction always tends to suppress 
shot noise. 

It is also worth pointing out that our results ([25]) and 
(|27|) can be rewritten in terms of the following equations 



R 

Ri 
R 



R2 
R 



3/3-1 



(32) 



RA _ 157 -f5(3;g- 1)^-2 
R J ^ 18 



Two real solutions of these equations (corresponding to 
i?i ^ R2) can be easily found numerically. Thus, infor- 
mation about the first three current cumulants is suffi- 
cient to determine all three resistances Rd and R2 
in the system under consideration. 

Finally, we would like to note that our analysis also 
allows to derive the fourth and even higher current cu- 
mulants for our composite conductor. For this purpose 
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it is necessary to establish higher order in x terms in the 
expression for the generating function iS{x)- This cal- 
culation, although quite tedious, can be performed in a 
straightforward manner along the same lines as it was 
demonstrated here. Yet another promising extension of 
our formalism could be to apply it to hybrid normal- 
superconducting structures. In this case modifications 
simply amount to including superconductivity into the 
Usadel equations in the standard manner. The corre- 
sponding analysis, however, is beyond the frames of the 
present paper. 
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Appendix A 

Let us present some details of our calculation. The ma- 
trix ge{0) satisfying the normalization condition g^{0) = 
1 can be written up to the third order in x as 



-fqs 



1 2/ 

-1 



93 



-fP3 P2 
P3 fP3 



' P2P3 ) 



fq3 



92 



- P2P3 ) 



-fr-3 -B r2 
r3 frz + B 



(Al) 



where 



B = fpiqi + - {fp2pl + fpl + P392 + ^293) ■ (A2) 



Resolving the self-consistency equations and (|24p . 
for / and we obtain 



^^^i.(A3) 



The expressions (IA3|) account for the first order in x 
term in the generating function iS{x)- The contributions 
iIxRi,2/cR come from the interface terms dTTI) . while the 
term UxRd/cR emerges from Eq. ([T0| describing a dif- 
fusive element conductor. Their sum yields the first term 
in Eq. i.e. just the Ohm's law for our composite 

conductor. 

Next let us reconstruct the second order in x contribu- 
tion to the generating function iS{x)- For that purpose 
we need to determine the parameters p2 and 93. They 
are 



P2(e) 



2iLpRi 
3R 



N^{e-eV) 2 - — i 



3i?i Rl + Rl 



N'it) 1 



iipRi 

IT' 



93(e) 



2R 



N{e -eV) + - iN{e) ~ N{e - eV)) 1 + 



Rf + Rl 
i?3 



(A4) 



2 {Rl + Rl 
i?3 



The zero-frequency current noise follows from the second 
term in Eq. (^5]). 

Finally, in order to find the third current cumulant 
it is necessary to find the parameters (72 and r^. The 
corresponding calculation is straightforward but yields 
rather lengthy expressions which we do not present here. 
One encounters the following integrals 



dx (tanh x — tanh(a; — a)) tanh'^ x 



2a 4a , 

= \ F(2a), 

3 3 ^ ^' 



dx (tanh x — tanh(a: — a)) tanh x tanh(a; — a) 

2a 8a , ^ , , ^ 

- - -F{2a), (A5) 



which contain the function F{v) (|28p . Collecting all 
terms one eventually arrives at the third order in x con- 
tribution to the generating function p5|) - (|27l) . 
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